Optimal Measurement on Noisy Quantum Systems 
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We identify the optimal measurement for obtaining information about the original quantum state 
after the state to be measured has undergone partial decoherence due to noise. We quantify the 
information that can be obtained by the measurement in terms of the Fisher information and find 
its value for the optimal measurement. We apply our results to a quantum control scheme based on 
a spin-boson model. 
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The most serious obstacle against realizing quantum 
computers and networks is decoherence that acts as a 
noise and causes information loss. Decoherence occurs 
when a quantum system interacts with its environment, 
and it is unavoidable in almost all quantum systems. 
Therefore, one of the central problems in quantum in- 
formation science concerns the optimal measurement to 
retrieve information about the original quantum state 
from the decohered one and the maximum information 
that can be obtained from the measurement. 

In this Letter, we identify an optimal quantum mea- 
surement that retrieves the maximum information about 
the expectation value of an observable of p from the par- 
tially decohered state. Here, p is an unknown quantum 
state and modeling of the noise is assumed to be given. 
The information content that wc use is the Fisher infor- 
mation [H , which has been widely used in estimation 
theory and is related to the precision of the estimation. 
For cases in which the unknown quantum state can be 
described by a single parameter, an optimal procedure 
to estimate this parameter has already been found [ij 
and used for phase estimation Q . In general, a quantum 
state is described by multiple parameters. The optimal 
estimation procedures for the multiparameter case have 
been discussed in several models of quantum systems Q 
and these are deeply connected with the uncertainty rela- 
tions of non-commutable operators [H . The main result 
of the present study is to identify the optimal measure- 
ment for a noisy quantum system (see also @). Here, by 
optimal, we imply that the Fisher information obtained 
by the measurement is maximal and that the precision 
of the estimation from the measurement outcomes is also 
maximal. While the aim of quantum error correction Q 
is to protect the unknown quantum state from interact- 
ing with the environment, our aim is to extract maximum 
information from the noisy quantum system. 

The crucial observation for obtaining our results is that 
the quantum state, observables, and quantum measure- 
ments are all described by a common set of generators 
of the Lie algebra. This fact greatly facilitates the anal- 
ysis carried out in the present study. The Fisher infor- 



mation describes the precision of the parameter estima- 
tion and it is defined through the parameterization of 
quantum states. We use a generalized Bloch vector Q 
as the parameter. Any quantum state of a finite iV- 
dimensional quantum system is expressed in terms of 
generators A = of the Lie algebra su{N). Let 

the generalized Bloch vector 8 G R^'~^ be defined as the 
coefficient vector of the expansion of /5 by A: 
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where / is the identity operator. Since p is unknown, 
6 is also unknown. The generators A satisfy aJ = Ai, 
Tr Ai = 0, and Tr[AiAj] = 25ij, and each A is character- 
ized by the structure constants fijk (completely antisym- 
metric tensor) and gijk (completely symmetric tensor) as 

[Ai,Aj] 2iJ2kfijk\k, {Ai,Aj} = ^Sijl + 2Y,k9i]k^k, 
where [ , ] and { , } denote the commutator and the 
anti-commutator, respectively. 

The quantum noise in a finite-dimensional quantum 
system can be described as an affinc map £ 0, S{p) = 
^iMipM] , where {Mi} are the Kraus operators that 
satisfy M^Mi = /. The Bloch vector 9 is also affine- 
mapped by £. By assuming that the dimension of the 
decohered state £{p) is the same as that of p, 

£{p) = ii+'^{Ae + c)-\, (2) 

where A is an (N'^ — 1) x (A^^ — 1) real matrix whose 
ij-element is iTr[AiiS(Aj)] and c £ whose ith ele- 

ment is ■^Tr[Ai£'(/)]. We assume that £ is injective [lo| : 
then, A has an inverse, which physically implies that £{p) 
is a partially (not completely) decohered state. The ob- 
servable X can also be expanded hy X as X ^ x^I + x ■ X, 
where G M and x G Then, the expectation 

value of X is calculated to be {X) ~ xo -\- x ■ 0. There- 
fore, estimating (X) is equivalent to estimating x -9, and 
our problem reduces to finding the measurement that 
maximizes the Fisher information about x ■ 9. 

We next introduce the Fisher information. Given n (3> 
1) independent and identically-distributed (i.i.d.) quan- 
tum states £{p), we perform the same POVM (Positive 
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Operator Valued Measure) measurement E = {Ei}i on 
each of them. The probabihty distribution of measure- 
ment outcomes is given by pi = Tr[£'(p)£'i]. In terms of 
Pi, the Fisher information about {X) obtained by E is 



defined as [11 j 



J{x;E) = [x-J{E)-^x]' 



(3) 



where J{E) is an {N'^ ^ 1) x {N'^ — 1) symmetric matrix 
called the Fisher information matrix, whose ij-element 
is defined as [JiE)]^j ^ Y^k^M^- Since J{E) 
has some zero eigenvalues, the inverse is defined on the 
support of J{E), which we denote as supp[J(£J)], and 
J{x; E) is defined as zero if a; ^ supp[J(£^)]. 

The Fisher information characterizes the precision of 
the estimation. The precision of the estimated value 
(estimator) X* of the unknown {X) can be measured 
by the variance oi X* . If the estimator X* satisfies 
the unbiasedness condition, that is, if the expectation 
value of X* for all possible outcomes equals {X), the 
variance Var(X*) satisfies the Cramer- Rao inequality: 
n Var(X*) > [J{x\ where n is the number of the 

samples that we measure. In general, the equality of the 
Cramer- Rao inequality is asymptotically satisfied for any 
POVM E by adopting the maximal-likelihood estimator 
as X* . Then, the estimation can be carried out most 
precisely with the measurement that maximizes J{x] E). 

The primary finding of our study is that the optimal 
measurement for obtaining the Fisher information about 
{X) is the projection measurement Py corresponding to 
the spectral decomposition of an observable Y that is the 
solution to the operator equation 



(4) 



where £^ {Y) = Y^t^lYM, is the adjoint map of £. 
Since Tt[£{p)Y] = Tr[pf the observable Y = yj + 

y-Ais adjoint mapped asf^(F) ~ {yo+y-c)I+{A^y)-X, 
where T denotes the transpose. Because we assume 
that A has an inverse, the solution to ([4]) is obtained 
as y = (xo - {[A^]-^x) ■ c) + [{A'^)-^x] ■ X. Although 
the Fisher information depends on the unknown quan- 
tum state p [l3l, the observable Y is independent of p. 
Therefore, Py is also independent of p, and the optimal 
procedure to estimate (X) is simply performing Py to 
the noisy system. We also find that the maximum Fisher 
information about {X) is given by 

Jix; Py) = {AY)-' EE {tt[£{p)Y'] Tr[f (p)f]'}"' . 

We can also use quantum state estimation strategies [l^ 
to estimate (X). However, these strategies provide un- 
necessary pieces of information about the system at the 
expense of decreasing the precision of the estimation of 
(X). Therefore, to estimate the expectation value of a 



single observable (X), performing Py is the best strat- 
egy- 

To prove these results, we first show that the Fisher 
information about (X) obtained by the projection mea- 
surement of Y is expressed as ([5|). Let P be a projection 
measurement. Because the elements of P = {Pi}iLi are 
Hermitian operators, they are expanded in terms of A as 
Pi = j^I+Vi-X, where Vi G M.^ For the completeness 

of the measurement, Vi must satisfy J2^=i = 0- When 
we measure £{p) with P, the probability distribution of 
the outcomes is given by Pi = ji + Vi ■ {AO + c) . Then, 
the Fisher information matrix J(P) is calculated to be 
J(P) = A^KA, where K ee Y.f^^P~^VivJ . To calculate 
the Fisher information about {X), we need to find the 
inverse of K. The support of K is the space spanned 
by {fij^j. The inverse of K for supp(i^) is given by 
if-i = (yT)-iQy-i^ where V is an {N' - 1) x ma- 
trix whose it\\ column vector is Vi, and Q is an x 
symmetric matrix whose ij-element is 5ijPi — PiPj . Be- 
cause V is not a square matrix, we denote as the 
generalized inverse matrix oi V . If we express the sin- 
gular value decomposition oi V as, V ^ Tli^id'Hi i the 
generalized inverse is defined as V^^ = s-^riiQj . 
We therefore obtain 



J{x- P) = [x- A-^K-^{A^)-^x\- 

= [v{v'^)-'QV-'v]-\ 



(6) 



for y = {A^)~'^x e supp(A:), and J{x;P) = for 
y ^ supp(7f). The condition y S supp(Ar) is equiva- 
lent to the condition that P is the projection measure- 
ment Py that corresponds to the spectral decomposition 
of an observable Y = y ■ X. By denoting the spectral 
decomposition of Y as Y ~ YiLi <^iPii it follows from 
the definition of V and the completeness conditions of 
P that the zth eigenvalue is equal to the ith element 
of V^^y G M^. Therefore, the Fisher information ob- 
tained from Py can be calculated to be the inverse of 
the variance of Y on £{p): 



J{x-Py 
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{AY)- 



(7) 

We next show that ([7]) gives the maximal Fisher in- 
formation. To show this, we use the quantum Fisher 
information [l3 | and the quantum Cramer-Rao inequal- 
ity [iB l. The quantum Fisher information matrix J'^ is 
independent of measurements, depends only on the mea- 
sured quantum state £{p), and gives an upper bound on 
the classical Fisher information matrix via the quantum 
Cramer- Rao inequality: 



J{E) < for aU E. 



(8) 



Therefore, the classical Fisher information J{x; E) is 
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bounded from above as 



J{x; E) < J'^ix), for all E and x, 



(9) 



where J^{x) = [x-{J^)^^x\^^ is the quantum Fisher in- 
formation about {X). Among several types of quantum 
Fisher information matrices that satisfy ([5]), we adopt 
the symmetric logarithmic derivative (SLD) Fisher in- 
formation matrix, which gives the tightest bound [HI on 



whose ij-element is defined as [J'^]ij = Tr 



dS(p) 



where Li is a Hermitian operator called the SLD op- 
erator. The SLD operator is given as the solution 
to ^£{p) = ^{E{p),Li}. Expanding the SLD op- 



erator 
b, = 



as 

-1 + 
£{0), where 
1, £(6) = AO 



Li = aj + bi ■ A, 

+ Gs<e)-£{e)£{eV) 



from 



m 

and a 



we obtain 
-6, • 



is a unit vector whose zth element is 
+ c, and Gffe) is a matrix whose ij- 
element is gijk£iS)k- From the definition of the 
SLD Fisher information matrix, its ij-element is cal- 



culated to be (Asi) ■ bj 



We 

-1 



thus obtain = 
^j^-. I ^£(^0) — £{6)£{9)^^ ^ A. Since we assume 
that A has an inverse, the SLD Fisher information about 
{X) is 



jQix) = [x ■ A-'i^I + Gs^o) - £{e)£{er){A^)-^x]-' 
= [Tv[£{p)Y']-Tr[£{p)YfY\ (10) 



Then, it follows from ([7]) that the projection measure- 
ment Py oi Y = y ■ \ satisfies the equality of © and 
that Py is the optimal measurement for obtaining the 
Fisher information about {X) from £{p). 

Since Ay and Py are invariant under transformation 

Y — !■ yo^+^" for any G we can choose the observable 

Y so as to satisfy £'^(Y) = X. Therefore, to estimate 
{X) from the decohered state £{p)^ the optimal method 
is to perform the projection measurement Py of Y that 
satisfies the operator equation £'^{Y) = X. 

As an illustrative application of our results, let us con- 
sider a situation in which a single qubit interacts with a 
heat bath of bosons [l7i . The total Hamiltonian Hq is 



where 6^ is the bosonic annihilation (creation) oper- 
ator of the heat bath. We assume that the state of the 
qubit and the bath is separable at t = and that the 
initial state of the qubit is p and that of the bath obeys 
the canonical distribution. The state of the qubit at t is 
calculated in the interaction picture to be 

£ip; t) ^ i(l + e-r°«)p + 1(1 - e-r°«)a,/5a„ (11) 

where TQ{t) increases monotonically from zero at t = 0: 
Toit) = 4 dcoDiuj) ^-"^r* coth . Here, D{u) 



is the spectral density function of the bath that we as- 
sume to take the form D{uj) = jUJ e~'^^'^'' , where Wc is 
the Debye cut-off frequency. Then, A{t) and c{t) of 
this quantum operation (jlip are found to be A{t) = 
diag(e"^«(*),e-i'''(*), 1) and c{t) = 0. Therefore, A{t) 
has an inverse for i < +oo. At i = +cxd, the right 
singular vector corresponding to the non-zero singular 
value of A{t) becomes (0,0,1)"""; then, the Fisher in- 
formation about all but az vanishes. If we substitute 
X ~ sin 6'obs'5'£!; + cos6'obs'5'z, then the solution to (jH) is 



Yit) 



cos6'obsO'z, so that the opti- 



mal measurement for £{p; t) is the projection measure- 
ment P±{t) = \i± \{s\ne{t)ax +cosd{t)az), where d{t) 
satisfies tand{t) = e'""^*-' tan^obs- Thus, the measure- 
ment direction tilts toward the x-direction and eventu- 
ally converges to the a;-direction, as shown in Fig. [TJa) . 
Moreover, the information about X except for 6'obs = 
converges to zero; therefore, we cannot estimate any ob- 
servable except for ir^ at t = +oo (sec dashed curves on 

Fig.inb)). 

In the above example, the qubit is decohered and the 
information about the system decreases monotonically 
because of the effect of the noise caused by the interac- 
tion with the heat bath. It is known that the decoher- 
ence for the spin is suppressed by the spin echo tech- 
nique by applying a sequence of pulses [itI [i8| . In this 
case, however, A{t) is not diagonal, and the measure- 
ment direction is drastically changed. We consider the 
case in which the total Hamiltonian is -ff = Hq + Hr{{t), 
where Hri{t) describes the effect of the pulse irradia- 
tions HAt) = Efc^^''HO{cos[cJo(t - 4^')]<T. + 
sin[cJo(^ — t'p^)]^y\ with t^jf' = kAt + (fc — l)r and 
]/('=)(*) = V for < t < 4''' + T = ife and other- 
wise. Each pulse is applied from tp*^' to tk, and the time 
interval to the next pulse is At. Here, amplitude V and 
duration r are tuned to satisfy ^ = f ■ Figures [Ijc) 
and (d) show the change in the measurement direction 
n{t) = (sin 9 (t) cos (l){t), sin 9 (t) sin (j){t), cos 6 {t)) of the 
optimal measurement P±{t) = ^(^ ± n{t) ■ &) for ob- 
taining the information about X. The solid curves in 
Fig. [Ijb) show the maximum Fisher information about 
an observable. By applying the sequence of pulses, most 
of the lost information is recovered; thus, the decoherence 
is suppressed. 

Here, we compare our optimal method with the quan- 
tum state tomography strategy [31 . For the example 
described above, the Fisher information obtained by our 
optimal measurement is three times larger than that ob- 
tained by the measurement proposed in [l9| . This is be- 
cause the quantum state tomography strategy divides a 
given set of samples for use to determine three noncom- 
mutable observables, whereas our strategy use all of them 
to determine a single observable. 

In conclusion, we identified an optimal method for es- 
timating the expectation value {X) from a noisy quan- 
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FIG. 1: (Color) (a)Time evolution of 6{t) without pulse irra- 
diation for kaT = lOhujc- (b)Time evolution of the maximum 
Fisher information J about X with ^obs = 0.257r for p = iJ, 
At = 0.3lo-^, and r = 0.05At. The red (blue) solid curve 
shows the high (low) temperature case with ksT = lO^c 
{ksT = hojc) with the sequence of pulses, and the dashed 
curves shows the case without pulses, (c) and (d)Time evolu- 
tions of 0{t) and 0(t) of the optimal measurement when the se- 
quence of pulses is applied, where ksT — lO&Jc, At — 0.3a;~^, 
and r — 0.05 At. In (d), the time scale of pulse irradiation is 
magnified for clarity. 

turn system. The optimal measurement that maximizes 
the Fisher information is the projection measurement Py 
corresponding to the spectral decomposition of Y that 
satisfies £^{Y) — X. Wc also find that the maximum 
Fisher information obtained by the measurement is given 
by the inverse of the variance of Y for the decohered 
state. Although the Fisher information depends on the 
unknown quantum state, the optimal measurement that 
maximizes the Fisher information is independent of the 
unknown quantum state. Therefore, the optimal strat- 
egy for estimating (X) is to perform Py on the noisy 
quantum system. Our results are obtained under the 
assumptions that the quantum noise £ is injective and 
that the Hilbert space of the original state p and the de- 
cohered state £{p) have the same dimension. The non- 
injectivcness of £ corresponds to the case in which the 
quantum state is completely decohered by the noise, for 
example, at t ~ +oo in the previous example. When 
quantum states are transferred by or stored on other me- 
dia, we can envisage situations in which the dimensions 
of the Hilbert space of p and £{p) are not equal. There- 
fore, solving the problem in such situations is crucial for 
implementing quantum networks and memory. The full 



investigation of this study will be reported elsewhere. 
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